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Abstract 

We analyse the OPE contribution of gluon bound states in the double scaling limit of 
the hexagonal Wilson loop in planar AT = 4 super Yang-Mills theory. We provide a system¬ 
atic procedure for perturbatively resumming the contributions from single-particle bound 
states of gluons and expressing the result order by order in terms of two-variable polyloga¬ 
rithms. We also analyse certain contributions from two-particle gluon bound states and find 
that, after analytic continuation to the 2 —>■ 4 Mandelstam region and passing to multi- 
Regge kinematics (MRK), only the single-particle gluon bound states contribute. From 
this double-scaled version of MRK we are able to reconstruct the full hexagon remainder 
function in MRK up to five loops by invoking single-valuedness of the results. 
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1 Introduction 

Light-like polygonal Wilson loops have been the subject of much study recently, in particular 
because of their relation to scattering amplitudes in the planar limit of AT = 4 super 
Yang-Mills theory [1-9]. The conformal symmetry of the theory and its associated Ward 
identity [10] imply that the Wilson loops with four or five sides are essentially trivial, with 
no interesting dependence on the configuration of the contour. From six points onwards 
however there is a conformally invariant function of the loop contour which needs to be 
determined. Here we will focus on the bosonic hexagonal Wilson loop, which corresponds to 
the six-particle MHV amplitude under the amplitude/Wilson loop duality. The conformally 
invariant function of the loop is known as the ‘remainder function’ in this case and is a 
function of three cross-ratios. 

One approach to studying the kinenratical dependence of light-like Wilson loops is to 
make an ansatz based on the analytic behaviour of explicit results obtained for terms in 
the perturbative expansion such as [11-13]. This approach has yielded results up to four 
loops at six points [14-18] and up to three loops at seven points [19]. The ansatze are 
given in terms of iterated integrals over words formed from a specified set, or alphabet, 
of rational functions (known as letters). The alphabets appearing in the explicit two-loop 
expressions of [13] were observed [20] to be described by the M-coordinates of a class of 
cluster algebras [21, 22] associated to the Grassmannians G(4,n). Taking this observation 
as an assumption, together with basic analytic information on the locations of possible 
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cuts of the final expressions produces a rather restrictive ansatz for the relevant functions. 
Indeed, at three loops and seven points, these analytic assumptions were essentially used to 
replace the dynamical information of the theory, producing a unique result for the symbol 
of the heptagon remainder function [19]. 

A second very powerful approach to describing the Wilson loops is based on a type of 
operator product expansion applied to configurations of null Wilson lines [23-27]. In this 
approach the Wilson loops are given by an infinite sum over excitations of a light-like flux 
tube. This sum amounts to an expansion around a collinear limit, and the excitations may 
also be thought of as insertions of the fields of the theory (gluons, fermions and scalars) 
on the Wilson line segments that are becoming collinear. The number of excitations or 
insertions is equal to the order at which these appear in the near-collinear expansion. 
The spectrum of excitations was calculated exactly in [28] using techniques based on the 
integrable structure exhibited in the planar AT = 4 theory. Apart from the spectrum, a set of 
overlap functions, called pentagon transitions, are required to calculate the near-collinear 
expansion. In [29] an all-loop formula was proposed to describe the subset of pentagon 
transitions consisting of any number of gluons, including multi-particle bound states. 

It is interesting therefore to understand how the near-collinear OPE expansion can be 
resummed into the type of cluster polylogarithmic functions appearing in the perturbative 
expansions of the null Wilson loops. As a first step towards the resummation of the OPE 
series, we develop a systematic procedure for summing the contributions to the OPE of 
single-particle bound states of an arbitrary number of gluons in the weak coupling expansion 
of the MHV 6-point (or hexagon) remainder function. Since these states have a restricted 
dependence on the helicities of the gluons, this can really be thought of as a contribution to 
the ‘double scaling’ limit of the Wilson loop, where one of the three cross-ratios is taken to 
zero. In the perturbative regime the double scaling limit is entirely governed by the gluon 
insertions, i.e. to determine it entirely at each perturbative order one needs all multi-particle 
bound states of gluons but not the contributions of fermions and scalars [29]. 

Building on the previous works [30, 31], our procedure relies on the technology of 
nested sums [32], and in particular its nestedsums C++ library realisation [33]. Using 
these algorithms we are able to show that the OPE expansion for single-particle gluon 
bound states can always be resummed into two-variable polylogarithmic functions based on 
a five-letter alphabet. These functions may be described as +2 polylogarithms in the cluster 
algebra language [20] or polylogarithms on .Ado ,5 hi the language of [34] or as a subset of 
the two-dimensional polylogarithms of [35]. In fact we find that only a particular subset 
of these functions arises, which is consistent with the functions having restricted branch 
cuts, and indeed with the idea that they are particular limits of hexagon functions [16, 17] 
describing the full six-point remainder function. 

Perhaps more importantly, from the knowledge of the single-particle bound state con¬ 
tributions to the double scaling limit we can produce all but the power-suppressed terms of 
the MHV hexagon in multi-Regge kinematics corresponding to high energy gluon scattering 
in the 2 —> 4 Mandelstam region. The kinematic limit on the cross-ratios is formally the 
same as the soft limit, in which the remainder function vanishes. To obtain a non-trivial 
dependence in multi-Regge kinematics, it is necessary to analytically continue the remain- 
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der function to the 2 —> 4 Mandelstam region, achieved by going around a singularity where 
one of the cross-ratios vanishes before taking the limit. In [36] it was observed that this con¬ 
tinuation can be carried out order by order in the OPE expansion. Using this property we 
provide evidence that after analytically continuing and taking the limit, the single-particle 
bound states we have summed are the only excitations with non-vanishing contributions 
on a one-dimensional double scaling slice of the full two-dimensional space parameterising 
multi-Regge kinematics. In other words we find that, after analytic continuation, the con¬ 
tribution from multi-particle bound states of gluons in the double scaling limit is always 
suppressed when taking the multi-Regge limit. 

Under the well-justified assumption that the natural function space in these kinemat¬ 
ics is given by the class of single-valued polylogarithms (or SVHPLs) [37], we may then 
reconstruct the remainder function in the full parameter space, from its knowledge on the 
line obtained from the double-scaling limit. We have carried out this procedure explicitly 
up to five loops. We emphasise that in our approach the contribution of the scalars and 
fermion insertions in the OPE expansion is non-vanishing, even after analytic continuation 
and taking the Regge limit. Indeed one can already see this in the analysis of [36]. The 
matter contributions just vanish on the line we are considering, hence we are actually re¬ 
constructing their contribution from the gluon bound states by invoking single-valuedness. 

Our logic is very much in line with the analysis of [38] who also reconstructed the 
multi-Regge limit of the hexagon amplitudes in the 3 —> 3 Mandelstam region from the 
OPE contributions of just the single-particle gluon bound states. Indeed, starting with the 
finite-coupling integral expression for the same OPE contribution that we consider at weak 
coupling, and performing an additional analytic continuation in Mellin space, the authors of 
[38] arrived at finite-coupling expressions for the MHV hexagon in multi-Regge kinematics. 
As this second continuation is fundamentally non-perturbative, it is quite interesting that 
in our approach we can trade it for single-valuedness directly at the perturbative level. 

The paper is structured as follows. We begin in section 2 by describing the OPE 
expansion of the hexagon Wilson loop in the double scaling limit. Our focus is on the 
contribution of the single-particle gluon bound state contributions and their resummation 
order by order in perturbation theory into two-variable polylogarithms of a particular kind, 
which we perform in section 2.3. We then describe in 2.4 why the class of functions so 
obtained is consistent with the idea that the full hexagonal Wilson loop remainder function 
is expressed in terms of hexagon functions. 

In section 3.1 we describe the procedure we use to analytically continue the results of 
the resummation to the double scaling limit of the 2 —> 4 Mandelstam region and then take 
the limit to (double scaled) multi-Regge kinematics. Having obtained a particular limit of 
the result in multi-Regge kinematics we then describe in section 3.2 how we can complete 
our expressions to restore the full kinematical dependence in multi-Regge kinematics by 
demanding that the result is single-valued. The analysis of this section is also supplemented 
by appendix B, where several two-particle bound states are computed and shown not to 
contribute to MRK. 

Attached to the arXiv submission for this paper are hies containing our results for the 
resummation of the single-particle gluon bound states up to five loops, particular contribu- 
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Figure 1. Decomposition of the light-like hexagonal Wilson loop in to a top pentagon, bottom 
pentagon and intermediate square. The collinear limit is indicated by the arrows, and each term in 
the expansion around it may be mapped to an excitation of an integrable colour-electric flux tube, 
sourced by the two sides of W\j adjacent to the ones becoming collinear. 


tions from two-particle gluon bound states contributing to the double scaling limit and the 
new results at N 3 LLA and N 4 LLA for the completion of the five-loop remainder function 
in multi-Regge kinematics. 


2 The Hexagon Wilson Loop OPE 


2.1 Preliminaries 

An operator product expansion (OPE) for light-like Wilson loops was introduced in [23] 
and refined in many papers [24-27, 29, 39-42]. It describes the near-collinear regime of 
a particular ratio of light-like Wilson loops, denoted by W. Here we will focus on the 
hexagonal Wilson loop where the ratio takes the following form, 


w= WWa 

kktopHbot 


( 2 . 1 ) 


Here W is the hexagonal Wilson loop, while Wt op and W^ot correspond to the pentagonal 
Wilson loops and Wa to the square Wilson loop indicated in Fig. 1. 

The ratio W is finite and conformally invariant and is hence a function of the three 
available conformal cross-ratios, 


™2 „2 
_ X 46 rg 13 

u/1 2 2’ 

x 36 x 14 


„2 

_ x 15 x 24 
^2 2 2 ’ 

X 14 X 25 


U 3 = 
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~2 „,2 

x 25 x 36 


( 2 . 2 ) 


The OPE describes the ratio W in an expansion around the limit where the two adjacent 
edges X 56 and X6i, of the hexagon become collinear so that x f 5 vanishes. The limit is most 
conveniently parameterised by the variables {r, cf), a} which are related to the conformal 
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cross-ratios as follows 


1 e 2cr+T secliT 

Ul ~ 2 1 + e 2a + 2 e a ~ T cos </> + e- 2r ’ 

u 2 = -e _r sechr , (2-3) 

1 

3 1 + e 2a + 2 e a ~ T cos 4> + e _2r 

The limit r —> oo corresponds to the collinear limit. The quantity W becomes 1 in the 
r —> oo limit, up to exponentially suppressed corrections. The Wilson loop OPE gives a 
prediction for the form of the exponentially suppressed corrections, namely at l loops 

oo [rrx/2] l—l 

w® = £ e- mr E cos l< m - 2 2>W E • (2.4) 

m= 1 p =0 n =0 

Here [x] denotes the integer part of x. 

The basic idea behind the OPE is to express the bottom part of the Wilson loop as 
a coherent sum of excitations of the GKP string. These excitations then propagate to the 
top part of the Wilson loop where they are absorbed. Both the spectrum of excitations 
[28], which controls the propagation of states, and overlap functions [25], describing their 
production and absorption, can be studied at finite’t Hooft coupling using integrability. 
The OPE expansion then has the following schematic form 

W = P(0\i’)P(ip\0)e- E ^ T+ip ^ a+im ^ . (2.5) 

It is a sum over intermediate states t/j, weighted by the overlap functions for production and 
absorption (called ‘pentagon transitions’) P and a factor due the propagation from bottom 
to top involving the GKP energy (or ‘twist’) E^, momentum and helicity m^. For the 
purposes of this paper all the relevant quantities are available in the literature and we will 
describe them in greater detail in the following subsections. 

The ratio W is very simply related to the remainder function Rq via 

Re = log W - log Weds , (2.6) 


with 


log Weds = ^p{Li 2 (u 2 ) - Li 2 (1 - ui) - Li 2 (1 - u 3 ) + log 2 (1 - u 2 ) 

7T 2 

- log (td) log ( u 3 ) - log {ui uz) log (1 - u 2 ) + —} , (2.7) 

fa 

and where the overall coefficient T cusp in the latter formula is the cusp anomalous dimension, 
whose expansion in the coupling 


g 2 = 


(4tt) 


a 

2 


( 2 . 8 ) 
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is as follows 
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+ 8C3J 
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(2.9) 


In the next section we will describe the ‘double-scaling’ limit of W (and hence the 
remainder function Rq) which allows us to consider contributions to the sum over states 
(2.5) coming from gluons only. 


2.2 The gluon contributions and the double-scaling limit 

Apart from the collinear limit we described in the previous section, another kinematical 
limit which will relevant for our discussion is the so-called [24, 29] 

‘double scaling limit’: r, i(f> —> 00 , —r + i<f> fixed . (2.10) 


At the level of the OPE expansion (2.4) it amounts to the subset of contributions with 
p = 0 and cos mcf> —> e* m< ^, which are evidently the only ones surviving the limit. And as 
far as the cross-ratios are concerned, from (2.3) and (2.10) we deduce that 


1 


U2 


U 1 


U 1 = 


^ _|_ g— 2 cr _|_ g—<7 


^3 


^3 = 


2 _|_ g2 <7 _|_ gcr— t+20 ’ 


( 2 . 11 ) 


namely the double scaling limit describes a two-dimensional subspace of the kinematics 
with vanishing 112 , but general values for the remaining cross ratios. 

From the point of view of the OPE, this limit is interesting because a restricted set of 
relatively simpler flux tube excitations contribute to it. In particular, these are the gluon 
excitations with positive helicity, given that all fundamental scalar, fermion and gluon 
excitations have helicity (the charge conjugate to cj) hi (2.5)) 0,±l/2 and ±1, all of them 
have twist 1 at weak coupling, and bound states are formed between gluons with the same 
helicity, with their charges being just the sum of the charges of their constituents. It is 
important to note that the twists receive corrections at each order in the coupling while the 
helicities do not. It is therefore important that, in restricting to the gluon contributions 
only, we are first expanding perturbatively in the ’t Hooft coupling g 2 before taking the 
double scaling limit. That is, in the limit in question 1 


^6 —t Re = logWg + — log Wbds , (2-12) 

where Wbds is obtained from (2.7) after substituting (2.11), and the all-plus gluon con¬ 
tribution is a sum over an arbitrary number of effective particles N, each of which may 
consist of any number of bound gluons a*, i = 1,... ,N [29] , 


°o °o ai a N -1 * / \ ~ / \ 

yy _ 1 / dUi...aUN HaAUi) . ..il aN (u N ) 


N =0 ai=l <22=1 aj^ = l k^O 


nk ! 


(271-)^ n Pa ila M\uj) 

* 7 4? 


(2.13) 


1 Note that although Wbds has an unphysical branch cut at 112 = 1, as can be seen from (2.7), this is 

not visible in the double scaling limit, since 112 —> 0. 
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In the above formula rik counts the number of bound states made up from = k gluons, 

£a(«) = fia ( u )e- E ^ T+ip ^ u)a+ia ^ , (2.14) 

and the energy E 0 . momentum p a , and measure g a of the a-th gluon bound state, as well 
as its pentagon transition to the 6-th bound state P a \b, are reviewed in appendix A. 

As a final remark, notice that the weak coupling expansions (A.14)-(A.17) of all the 
different ingredients of (2.13)-(2.14) , imply that the IV-particle all-plus gluon state starts 
contributing at order 0(g 2N2 ). Namely one need only consider the states with N = 1 up 
to three loops, and IV < 2 up to 8 loops. 


2.3 Resumming all single-particle gluon bound states 

In this section, we will present a systematic procedure for perturbatively resumming the 
contribution of all single-particle gluon bound states to the hexagon Wilson loop/scattering 
amplitude, namely the N = 1 term in (2.13), which we may rewrite as 


Wi = 9 2l M l] 

l=i 



du 

2n 


H a {u)e~ Ea ( U ) T +*P“ {u)<r+ia(j> 


(2.15) 


This procedure is an extension of the method developed in [30, 31] for the evaluation of the 
individual a = 1 and a = 2 terms above. After describing its details, we will apply it in 
order to obtain explicit for expressions up to l = 5 loops, and also deduce the relevant class 
of functions describing this contribution to arbitrary loop order. The reader interested in 
the final result may jump directly to the discussion around eqs. (2.29)-(2.30) and (2.32). 

We start by expanding the integrand in (2.15) with respect to g, following a well-known 
procedure that we review towards the end of appendix A. This reveals that up to u- and 
a-independent factors (e.g. e~ T or powers of t and a) which we can factor out, all terms 
are of the general form 


J2(-e~ T+i *) a 

a=1 


du e 2iUCT T(f + i«)T(| - iu) n, V> (mi) (l + § ± ill) El i ± iu) 

27t r(a)( U + fY(u- f) r ' 


(2.16) 

for different powers of the denominators r, r ' and different products/powers and orders 
rrii^m'j of the polygamma functions ip, see eq. (A. 13) for their definition. 

From (2.16), we can immediately infer that at any loop order the poles of the integrand 
are located at u = ±i(a /2 + k), k = 0,1,2,.... After we restrict to a > 0, we can close the 
contour of integration on the u > 0 plane, and trade the integral for a sum over residues by 
virtue of Cauchy’s residue theorem. With the help of the recurrence and reflection relations 
of the gamma function, 


T(z + l) = zT(z), T(1 - z)T(z) = , (2.17) 

sin (7 tz) 
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as well as analogous equations for polygamma functions arising upon differentiation of 
(2.17), we arrive to sums over residues which always have the form 2 


» 

a= 1 


—e 


—T+i4 


oo 
<t>Y 

k=1 


T(a + k) Yli' l P {ni) ( a + k)Ylj^ { - n ' i \k +!) 


r(a)r(A: + l) 


(a + k) s k s ' 


(2.18) 


The next step is to reexpress the polygamma functions in terms of S- or Z-sums [32] via 3 , 


where 


tp(k + 1) = ^°\k + 1) = -7J5 + S{k- 1; 1) 
V' (m -i ) ( k + 1) = (-1 ) m {m - l)![Cm - S(k- m; 1)] , 


5(n; mi,..., rrij\x \,..., Xj) 
Z(n\ mi, ..., mj;x\, ..., Xj) 


E 


V 1 




•mi • ■ 


•j 


E 


•mi 


n>ii>i 2 >...>ij >0 1 


.rrij 

l- 


(2.19) 

( 2 . 20 ) 
( 2 . 21 ) 


and replace the products of S’- or Z-surns with the same outer summation index with linear 
combinations thereof, simply by nesting the independent summation ranges of each term 
in the product 4 . For example, 



n i\ 


n %2 


EE + EE- E 

Kh=l z 2 —1 i 2 =lii=l h=i 2=iy 
^ n ii — 1 n i 2 — 1 n N 

EE + EE+ E 

=1 *2 = 1 * 2=1 * 1=1 * 1=*2 = 1 2 


1 


•mi -m 2 
4 1 *2 

1 

•mi -m 2 > 

n '2 


( 2 . 22 ) 


where the first line yields S’-sums, and the second line Z-sums. 

For reasons that will become apparent very shortly, we will choose to replace ^ m \k) 
by S- and i// m )(a + k ) by Z-sums respectively. After shifting the summation variable 
a —> j = a + k, partial fractioning in k and expanding, (2.18) splits into terms that look 
like 


OO 


E 


_ p —T+i(j>—a 


(-e 


y 


j -1 


o'ni 


^0'- l;n 2 ,...;l,...,l)^ 


fc=i 


j - 1 
k 



k 

-S(k-ri 2 , 


-; i, - - -, i), 

(2.23) 


where we have combined the gamma functions of (2.18) into a binomial coefficient. 

Quite remarkably, the sum in /c can be done for any collection of n,, with the help of 
algorithm C of [32], Let us give a simple example where n\ = 1 and the S'-surn is absent, 

2 The k = 0 residue is treated separately, as it is the only case where the denominators have poles, and 
yields simple sums which may be evaluated exactly as in [30]. 

3 In particular, S{k\m\ 1) = Z(k\m\ 1) are the generalised harmonic numbers, the Riemann zeta 
function, and 7 e = —V>(1) =— 0.577 the Euler-Mascheroni constant. 

4 This is nothing but the quasi-shuffle algebra property of these objects. 
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so as to convey the basic idea of the algorithm, which is to appropriately manipulate the 
expression so that the sum in k can be done by means of the binomial theorem, 

£("V = (1+*)". (2*24) 

k =1 ' ' 


For the example in question, this entails getting rid of the denominator, roughly speaking 
by rewriting the summand as the integral of its derivative with respect to x, 


n / \ l. 

n\ x 


E 

k =1 


k k 


SCOjCWS^- 1 ]- (-, 


By changing the integration variable to y = 1 +x, the integrand becomes equal to the sum 
of the first n terms of a geometric series, so that 


n 


E 


/ n\ x k 
\k)~k 



dy 


i -y n 
i -y 


n 



1+x 

1 


n 


E 


(1 + .x) fe - 1 

k 


S(n; 1; l+x)-S(n; 1; 1), 


where in the last line we used the definition (2.20). By recursively going through the exact 
same steps, we may similarly obtain 



dxm-i f x2 dxi _1_ 

Xm-l ’"Jo X 1 1 • ^ h 


= S(n; 1 ,..., 1; 1... 1,1 + x) - S(n ; 1,..., 1; 1... 1,1). 


1 (1 + x) lm - 1 

^m— 1 'I'm 


(2.26) 


In practice, the simplest way do any sum in k as in eq. (2.23), is to exploit an al¬ 
ready existing C++ implementation of the algorithm, as part of the nestedsums library 
[33] within the GiNaC symbolic computation framework [43] 5 . The relevant command is 
transcendental_sum_type_C, and we have built an interface that calls it directly from 
Mathematica, which we use for the remaining manipulations. Combining it also with the 
Ssum_to_Zsum command, the evaluation returns Z-sums with outer summation index j — 1. 

Then, after employing quasi-shuffle algebra relations generalising (2.22) in order to 
eliminate the resulting products of Z-surns in favour of their linear combinations, we may 
immediately evaluate the remaining sum in j in (2.23), in terms of multiple polylogarithms, 

OO A 

Li mi ,..., m .(xi, ...,Xj) = 22 - l;m 2 ,.. ■,m j ;x 2 , •. .,Xj), 

2=1 

see appendix C for more details on their definition and properties. 

The final result is most conveniently expressed in terms of the variables 

I = = “1 1+e r-<*-, = 

U\ 1 — U\ — Uo 

5 Alternatively, the XSummer package [44] for the FORM symbolic manipulation system [45] 
the same functionality. 


(2.28) 
offers exactly 


(2.27) 
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and for example at one loop we have (we denote (x) = Li1,..., 1) the 

subclass of harmonic polylogarithms [46]), 

wj 1) = —Li ljl (x, y) + 2aHi(x) + H 2 {x ), (2.29) 

and at two loops, 

Wf ] =t {a [—4Lii i i(x, y) - 4 H 2 (x)\ + 4Li 1> i j i(x, y, 1) - 4 H 3 (x)} + 4a 2 Hi A (x) 

+ cr — 4Lii ) i i i(x, y, 1) + ^(- 2)n 2 Hi(x ) -4iL 3 (x)^ (2.30) 

+ ^tt 2 H 1} i(x) - 2H lj3 (x) - 2H 2 fl{x) - 2H 3)l (x) - 2H 1 ^ 2 (x) - 2H 1)2 ,i(x) 

- 2iL 2 ,i,i(x) + 2Lii >3 (x, y) + 2Li 2 , 2 (x, y) + 2Li 3 p(x, y) + 2Lii i i j2 (x, 1, y) - GH^x) 
+ 2Lii )2 ,i(x, 1, y) + 2Li 2) i ! i(x, 1, y ) + 2Lii i i : i ) i(x, 1,1, y) + 4Lii.i ) i i i(x, y, 1,1). 

As we mentioned at the beginning of this section, we have computed up to l = 5 loops, 
and the remaining results may be found in the accompanying file allboundstates .m so as 
to avoid clutter. 

We stress that the procedure we have described can be applied in principle at any loop 
order, and in fact we can make a precise statement about the particular class of multiple 
polylogarithms in which lies: From (2.23) it is evident that the leftmost argument of 
the MPLs (2.27) will always be x. Furthermore, by inspecting the steps of algorithm C in 
[32], we observe that it can only generate Z-surns with arguments y,l/y and 1, such that 
y and 1/y always appear in alternating order with y leftmost, if one removes all arguments 
equal to 1. Combining these two requirements, we infer that only MPLs of the form 

Lim 1 ,m.2,... (a:, 1, — , 1, y, 1,..., 1, 1/y, 1,..., 1, y ,...) (2.31) 

can appear in our results. 

The basis of MPLs (2.31) we have obtained may be expressed more transparently in 
terms of G- functions, whose definitions we review in appendix C. More specifically, the 
equivalent Li- and G'-function representations of MPLs are related by (C.6), and it also 
proves advantageous to rescale all resulting G'-functions so that the rightmost argument 
becomes x as in (2.28), by virtue of the identity (C.7). 

In this manner, we finally arrive at the following important conclusion: Apart from 
explicit factors of a = —^ log(x(l — y)) and r ~ ^ logn 2 (in the double scaling limit), the 

OPE contribution of all single-particle bound states is expressed in terms of 

G(a\, ..., a n \x) , at E {0,1, 1/y} with a n ^ 0 , and aj ^ 1/y if aq = ... = (ij-\ = 0 , 

(2.32) 

at any loop order l. This basis of MPLs has dimension at weight n 6 , and in fact it turns 
out to be part of the two-dimensional harmonic polylogarithms (2dHPL) of Gehrmann and 

6 We have a total of 2 • 3 n_1 G-functions with a n ^ 0, minus 1 + X )^=3 2 • 3 fc-3 = 3 n_1 of them with 
a\ = ... = aj -1 = 0 and a,j = 1/y. 


10 



Remiddi, which were first introduced in the computation of four-point functions of three 
on-shell and one off-shell leg at two loops [35]. 

In the next section, we will explore whether these results are consistent with the fun¬ 
damental assumption of the hexagon amplitude bootstrap [14-18] for the class of functions 
describing six-particle scattering in general kinematics. 

Calculating the contribution of multi-particle bound states to the double scaling limit 
is considerably harder than computing the single-particle bound state contributions. We 
will leave the resummation over such states to future work. However we can say something 
about the OPE expansion of certain contributions and this is presented in appendix B. 

2.4 Hexagon functions in the double-scaling limit 

An important finding of the previous section, is that a particular OPE contribution (2.15) to 
the hexagon remainder function Rq(u\, U 2 , u 3 ) (2.6) is always expressed in the basis (2.32) 
at any order in the weak coupling expansion. More specifically, this contribution is the only 
one that survives in the double scaling limit (2.11) up to 3 loops, and is supplemented by 
the N = 2 particle gluonic contributions up to 8 loops, as we noted at the end of section 2.1. 
Given that Rq is believed to lie in the space of hexagon functions for general kinematics, 
here we would like to explain why our results are indeed consistent with the double scaling 
limit U 2 —> 0 of hexagon functions. 

First of all we note that hexagon functions are just a special subset of all iterated 
integrals based on the nine-letter alphabet {iij, 1 — m, y*} for i = 1, 2, 3 described in [14], In 
other words they are a subset of polylogarithms on .Mo, 6 > the moduli space of six points on 
a Riemann sphere. It is simple to see that in the double scaling limit 112 —> 0 these reduce to 
iterated integrals on the five-letter alphabet' {ui, u 3 ,1 — m, 1 — 113,1 — u\ — u 3 }, or in terms 
of the variables x and y defined above, on the alphabet {x, y, 1 — x, 1 — y, 1 — xy}. Thus 
the double scaling limit will give poly logarithms on Afo, 5 - A basis for such poly logarithms 
is given in terms of products of iterated integrals of the form G{m y ; y)G(m x \ x) where the 
weight vector m y is made of elements drawn from the set { 0 , 1 } and the weight vector m x 
is drawn from the set {0,1, ^}. Here we have chosen a specific contour of integration over 
integrable words in the five letters such that one performs the integration along the y-axis 
first and then in the x-direction. We could of course have chosen to do it the other way 
round, yielding a basis with x and y exchanged. 

Hexagon functions, as defined in [16], are a special subset of polylogarithms obeying 
restrictions on their branch cuts, namely that the locations of the branch points in the 
Euclidean region occur only at the boundary rij = 0 or Ui = 00 . At the level of the symbols 
of hexagon functions this is reflected in the fact that the initial entries are only the Uj and 
not any of the other nine letters. When taking the double scaling limit, the allowed branch 
point locations take the form 

~ 1 - x{l -y) 

u\ u\ = --, u 3 u 3 = — -. (2.33) 

1 — xy 1 — xy 

'The fact that hexagon functions reduce to this alphabet in the double scaling limit, and thus approach 
a subset of 2dHPLs, was noted independently in [47]. 
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In other words there should be no branch cuts around y = 0 or x = 1 at all and the branch 
cuts around y = 1 must match those around x = 0. 

To make the branch cuts around x = 0 explicit we may use the shuffle relations to 
extract any trailing zeros from the weight vectors of the G-functions with argument x. In 
other words we reexpress the G-functions with argument x as linear combinations of G'- 
functions which have no trailing zeros in their weight vectors and explicit factors of log x. 
We may then replace the explicit logarithms of x with logarithms of the product x(l — y) 
at the cost of redefining the G-functions with argument y for each element of the basis. 
However once we have done so we are no longer allowed any further branch cuts at y = 0 
or at y = 1. Thus we must have no further factors of the form G(m. y ;y) at all. 

Thus we find that a basis compatible with the constraints on the branch points at 
y = 0,1 is given by all products of the form 

log p (x(l - y))G(m x ; x) , (2.34) 

where there are no trailing zeros in the weight vector m x . 

Of course finding a function in the above basis does not necessarily satisfy the remaining 
constraint, namely that there should be no branch point at x = 1. Only special linear 
combinations of functions in the above basis will also satisfy this additional constraint. We 
have checked that the combinations appearing in the results derived above do indeed obey 
the property of having no branch points at x = 1. 

Finally let us mention that, as described in (2.32), we never have a weight vector in 
the G-functions which begins with a string of zeros (of any length) followed by a We 
have verified, by analysing the known four-loop remainder function [17], that this property 
holds also for the full in the double scaling limit. In other words it holds not only 
for the single-particle gluon bound state contributions but also for the two-particle ones. 
This property appears to be an additional constraint on the form of the remainder in the 
double scaling limit which, at least in the case of the single particle states, we can observe 
as coming from the structure of the OPE. 

3 From Collinear to Multi-Regge Kinematics 

3.1 The soft limit, analytic continuation and multi-Regge kinematics 

Having obtained expressions for the single-particle gluon bound state contribution to the 
double scaling kinematics we may now ask about relevant physical limits of our expressions. 
Since our aim will be to obtain the six-particle scattering amplitude in multi-Regge kine¬ 
matics we will briefly review the relevant kinematic limits. We recall that the soft limit and 
the multi-Regge limit are formally identical in that we send one cross-ratio (113 say) to one 
and the other two ( u\ and U 2 ) to zero so that the ratios 

u\ _ _1_ u 2 _ ww* 

1 — U 3 (1 + w)(l + w*) ’ 1 — 113 (1 + ir)(l + w*) 

are fixed. In the formulas above we have introduced the variables w and w* which parametrise 
the possible remaining dependence in the limit. In the soft limit we expect the remainder 
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function to simply vanish. However, after we analytically continue the amplitude to the 
2 —> 4 Mandelstam region (obtained by continuing around the singularity at U 3 = 0 via 
U 3 —>• U 3 e~ 2m ), the above limit is non-trivial and we expect a series of divergent logarithms 
whose coefficients are functions of the remaining variables w and w*, 

l -1 

—» (2vri) ^log r (l - u 3 )[g( l \w,w*) + (2iri)h^(w,w*)\ . (3.2) 

r=0 

The limit (3.2) has been studied in many papers (see for example [48-58]). 

In our preceding discussion from section 2 we were considering the double-scaling limit 
where one of the cross-ratios, 112 was taken to zero. This means that starting from the 
double scaling limit we actually only have access to multi-Regge kinematics in the regime 
where w* —> 0 with w fixed (or, switching the helicities of the gluon bound states, to 
the regime w —> 0 with w* fixed). We refer to this regime as double-scaled multi-Regge 
kinematics. 

We may relate the variables w and w* to the variables o, r and cj) via, 


w = re 




w = re 


r = e 


(3.3) 


Note also that the variable x introduced to parametrise the double scaling limit is related 
to w via x = — w and that therefore w* = —xe~ 21 ^ and hence terms power suppressed in 
the double scaling limit correspond to terms power suppressed in w*. 

Now let us turn to the properties of the expressions we have obtained in the various 
physical limits discussed above. An obvious property, manifest from the form of the OPE 
expansion is that all vanish in the strict collinear limit r —> 00 . Likewise, the expres¬ 
sions we have obtained vanish in the soft limit x —> 0 with y fixed. All of the G-functions 
of argument x are power suppressed as x —> 0 and there is one such G in each term of our 
results. 

The other soft limit y —> —00 with x fixed is related by the symmetry u\ o U 3 to 
the first. Our expressions also vanish in this limit. To see this however requires taking the 
limit y —> —00 on each of the G-functions appearing in our expressions. It is instructive to 
do this explicitly as we will need to perform exactly the same operation when taking the 
Regge limit. 

In order to carry out the limit y —> —00 on our expressions it is first useful to express 
them in such a way that no G-functions are divergent in the limit. To arrange this we use 
the shuffle relations to unshuffle any zeros or ^ appearing at the end of the weight vector 
so that each G-function either has a weight vector ending in a one or there are no ones on 
the weight vector at all. In the latter case we may rescale the argument by y so that these 
G-functions have weight vectors with entries 0 and 1 and argument xy (i.e. they are HPLs 
of argument xy). For example we write 

Gi 1 ’ 0, l', x) = -G(l; xy)G(0, 1; x) + G( 1; x)G(0,1; xy) + G(±, 0,1; x). (3.4) 

The limit y —»• —00 may now be taken straightforwardly using e.g. the HPL mathematica 
package [59] to give the limits of the G-functions with arguments xy as y —> — 00 . In this 
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way we find that all our resummed single-particle bound state contributions to the double 
scaling limit vanish in the soft limit y —> —oo. 

The very same limit is required to analyse the nrulti-Regge limit of the double scaling 
limit of W. However to obtain a non-vanishing contribution we must first analytically 
continue to the Mandelstam region by passing round the branch but starting at 113 = 0. 
This is easily achieved since all of the cuts at 113 = 0 are manifest in the form of our result 
since they appear as explicit logarithms of x(l — y) = 113/ui. We simply need to continue 
these so that logx(l — y) —> logx(l — y) — 2ni. In so doing we obtain the result for W in 
the Mandelstam region in the double scaling limit. 

Having analytically continued our results we may then go to multi-Regge kinematics 
by taking exactly the same soft limit as above, namely y — > —00 with x fixed. This time 
of course we obtain a non-vanishing result, dependent on x = —w. Finally, we express our 
results in a way that will be convenient for us to complete our expression in double scaled 
multi-Regge kinematics to the full multi-Regge kinematics. In order to do this we note that 
in the double scaling limit we have 

T ~~\ lo & M 2 = lo g(! - U 3 ) - ^ log (ww*) + ^ log(l + w) (3.5) 

and 

cr = -r - - log (ww*) = - log(l - u 3 ) - - log(l + w). (3.6) 

We may use these relations to rewrite all contributions in terms of the variables w, w* and 
(1 — 113) which parametrise multi-Regge kinematics. In fact we also rewrite log (ww*) = 
log(— w) + log(— w*) and neglect the divergent logarithm log(— w*). This is because these 
divergent logarithms are tied to the finite parts by the completion to single-valued polylog¬ 
arithms that we describe in the next section. 

The above discussion has focussed on the contribution of the single-particle gluon bound 
states to double-scaled multi-Regge kinematics. We have also analysed certain contributions 
from two-particle bound states, as described in appendix B. We find that after analytic 
continuation and taking the limit these contributions are all power suppressed. Hence, 
based on this preliminary analysis it seems that the single-particle states give the only 
non-vanishing contribution to double-scaled multi-Regge kinematics. 

3.2 Completion to full multi-Regge kinematics from single-valuedness 

Here we would like to explain why deriving the result in multi-Regge kinematics from the 
double scaling limit is sufficient to reconstruct the full result in multi-Regge kinematics. 
We know that in the multi-Regge limit of the remainder function in the 2 —> 4 Mandelstam 
region we will obtain an expansion in powers of the divergent logarithm which we make 
take to be log(l — 113 ) 8 , 

l -1 

Rf —» (2ni)^log r (l-u 3 )[g( l \w,w*) + (2iri)hl l \w,w*)]. (3.7) 

r=0 

8 Note that in this section we have changed the expansion parameter g 2 a = 2 g 2 to reflect most 
commonly adopted conventions in the subject. 
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The coefficients of the divergent logarithms are separated into an imaginary part gP and 
a real part hP . For each l and r both gP and lip are single-valued polylogarithms (or 
SVHPLs) [37], 

The expressions we have obtained from the analytic continuation of the double-scaling 
limit allow us to obtain the limit of each gP and hP as w* —> 0 with w held fixed 9 . We 
also discard any divergent logarithms of the form logu>* as we take this limit, keeping only 
the finite term. From this data we may reconstruct the full dependence of all contributions 
to the real and imaginary parts in multi-Regge kinematics by invoking single-valuedness of 
the result. 

In order to explain this point we briefly recall the construction of single-valued polylog¬ 
arithms given in [60]. We begin with Knizhnik-Zamolodchikov equation in a single complex 
variable z, 

Tz l{z) = P + tP) l{z) - (3 ' 8) 

Here eo and e\ are two free non-commuting generators. We take the solution Lq of (3.8) 
normalised so that as z -> 0 we have Lq(z) ~ Lfj n (z)z eo where Lfj n (z) is analytic around 
z = 0 and LQ n (0) = 1. The solution can be represented as a sum over all words in eo and 
e\ with coefficients which are harmonic polylogarithms (i.e. regularised iterated integrals 
over dlogz and dlog(l — z) away from z = 0), 

L 0 {z) = YmH m {z). (3.9) 

m 

Here the sum is over all words m in eo and ei. 

A second solution of the same equation may be taken to be 

Lip) = To(l — ~)|eo—»—ei, e\— >— eo ■ (3.10) 

It is normalised so that L\(z) ~ Lf n (z)(l — z)~ ei as z —> 1 with Lf n (z) analytic at z = 1 
and Lf n (l) = 1. 

The two solutions are related by parallel transport via a constant series $ (the Drinfeld 
associator), 

L 0 (z) = L^z)^ . (3.11) 

Here <h is a series in eo and ei with coefficients which are (shuffle-regularised) multiple zeta 
values, 

d> = ^mCm(m). (3.12) 

m 

Explicitly expanded as a series in eo and ei we have 

$ = 1 + [eo, e\]C ,2 + ([eo, [eo, ei]] — [ei, [eo, ei])C3 + • • • (3.13) 

We may now consider the analytic continuation of the solution around the singular 
points at z = 0 and z = 1. Under an analytic continuation around z = 0 (via z —> ze 27Ti ), 
the monodromy of the solution Lq is explicit, given the asymptotics as z —> 0, 

MqLq(z) = Lo(z)e 2nie ° . (3.14) 

9 We may treat w and w * as independent complex variables to perform this limit. 
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The monodromy of Lq around z = 1 may be obtained by transporting to z = 1, where the 
monodromy of L\ is explicit, and then back again, 

MiL 0 {z) = L Q {z)<S>- x e- 2me '<b . (3.15) 

Now, to construct single-valued polylogarithms in a complex variable z, we consider a second 
‘primed’ alphabet e' 0 and e'i and form the following series in all four letters eo, ei, e^e^, 

£(z,z) = L 0 (z)Lq(z) (3.16) 

The second factor is built on the primed alphabet while the symbol ‘~’ means that the 
words in e' 0 and e\ are reversed with respect to the ones in eo and ei appearing in Lq(z). 

Now the series C(z,z) will be single-valued if it has no monodromy around z = 0 or 
z = 1. The monodromy around z = 0 is given by 

MqC{z, z) = L 0 (z)e 2 ^ ieo e-‘ 2n <L , 0 {z). (3.17) 

The series will be unchanged if 

e o = e o ■ (3-18) 

Similarly the monodromy around z = 1 is given by 

MiC(z,z) = L Q (z)$- 1 e-^ 1 $&#^(&)- 1 L' 0 (z), (3.19) 

where d* 7 is built on the letters e' 0 and e\. The series is unchanged if 

$-i e - 27rie i$l>'e2™ e i (i')- 1 = 1, (3.20) 

which, together with (3.18), defines the e' 0 , e\ alphabet in terms of eo and ei. The individual 
coefficients of words in eo and ei in the series C(z, z) are then necessarily single-valued 
polylogarithms C m (z,z), 

£(z, z) = m£ m {z, z). (3.21) 

m 

In other words, at a given weight, there are exactly as many single-valued poylogarithms 
in z and z as there are harmonic polylogarithms in a single variable z built on dlog .2 and 
dlog(l — z). 

Given the asymptotics of Lq defined above, the limit of each £ m (z,z) as z —> 0 (now 
treating z and z as independent complex variables), is a series in divergent logarithms 
log z with the coefficient of the finite term being simply H m (z). In other words there is a 
unique completion of a given harmonic polylogarithm H m (z) to a single-valued polyloga¬ 
rithm C m (z, z) such that the finite term in the limit 7 —> 0 is H m (z). This fact was already 
noticed and found to be very useful in [61] in reconstructing the full form of the three- 
loop ‘Easy’ integral appearing in the correlation function of four stress-tensor multiplets in 
J\T = 4 super Yang-Mills theory. 

Thus from the above discussion it is clear that to complete our expressions given in 
terms of HPLs with argument (—w) into single-valued expressions with the correct limit, 
we simply replace each instance of H with C. In this way we obtain the full expressions 
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5-loop N 3 LLA contribution 



Figure 2. Log-linear plot of the imaginary part of the 5-loop hexagon remainder function in multi- 
Regge kinematics at (Next-to) 3 -Leading Logarithmic approximation, g[°\ on the line w = w*. 


5-loop N 4 LLA contribution 



( 5 ) 

Figure 3. Imaginary part g 0 ' of the 5-loop hexagon remainder function in multi-Regge kinematics 
at (Next-to) 4 -Leading Logarithmic approximation on the line w = w*. 

for fjr ^ and hr ^ dependent on w and w*. We have verified that our results up to five loops 
reproduce the known expressions derived in [17, 37]. In addition we have obtained the 
remaining five-loop terms g^’' and g q 5 \ We also obtain but the real parts may be 
simply related to the imaginary ones so this does not constitute independent data. We give 
plots of the new data g | ' and g 0 ' along the diagonal line w = w* in Fig. 2 and Fig. 3. 
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3.3 Comparison with BFKL and the approach of Basso, Caron-Huot and Sever 

From the reconstructed expressions derived in the previous section one can compare to the 
BFKL formula for the amplitude in multi-Regge kinematics, 


e Jfe-H** = COS7rWa6 + ^ (-1 r 


w\ f r0 ° 


ur 


L 


du 


-oo V 2 + « 2 / 4 


W | ^reg(^, n) 


x exp 


-u(v, n) ( iri + log(l - u 3 ) + * log —^ . 

2 1 + w r 


(3.22) 


The formula expresses the amplitude in multi-Regge kinematics as a Fourier-Mellin trans¬ 
form of a factorised expression involving the BFKL eigenvalue uj(is,n ) and impact factor 
<J> reg (i/, n) which encode all the kinematical dependence of the amplitude. The other quan¬ 
tities in the above formula are given by 


u a b = ~T cusp \og(ww *), 

O 



ww* 

(1 + w) 2 (l + w*) 2 ' 


(3.23) 


In principle, from our reconstructed expressions, we can then find $ reg (i/, n) and w(z/, n) 
so that they are consistent with the perturbative expansion of the amplitude in multi-Regge 
kinematics that we reconstructed from the double-scaling limit after analytic continuation. 
However, since an all-order form for these quantities was obtained in [38] following a dif¬ 
ferent (though similar) logic, we can simply compare our expressions for the amplitude in 
multi-Regge kinematics with those obtained from their formula. Doing so we find perfect 
agreement up to five loops 10 . 

Notice that the reconstruction argument of this section essentially amounts to the claim 
that the full BFKL expression can be reconstructed just from the knowing the residue at 
the first pole on the positive u axis, assuming that the perturbative expression is given 
in terms of single-valued polylogarithms in w and w*. Certainly, under this assumption, 
the simplest way to evaluate any given expression perturbatively is simply to calculate the 
integral above around the first pole and resum the result as a single-variable function of 
w. Then one may complete this boundary information to a single-valued polylogarithm as 
described above. 

It is very interesting that both here and in [38] the crucial ingredients were the single¬ 
particle bound states of gluons appearing in the OPE expansion. We stress that the other 
states do contribute in multi-Regge kinematics. It is just that we are able to ignore the 
contributions from the fermions and scalars from the beginning by going to the double 
scaling limit and then, by extension of our observations, assume that the multi-particle 
gluon bound states drop out when going to double-scaled multi-Regge kinematics. The 
contributions of the missing states can then all be reconstructed by appealing to single¬ 
valuedness. Note that the issue of identifying the n = 0 term in the sum (3.22) does not 
arise in our approach; it follows along with all the others when completing the double-scaled 
MRK expressions to the full MRK ones. If we were able to fully justify the fact that the 
multi-particle gluon bound states drop out from double-scaled MRK we would have proven 

10 Note that in the weak coupling expansions derived in [38], one needs to take care that r'BCS = 2 v. 
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that the full expression can be reconstructed from the single-particle gluon bound states 
alone. As the authors of [38] stress, very similar assumptions have to be made also in their 
approach. It would be very interesting indeed if some combination of the two arguments 
could be used to fully prove the BFKL formula for the hexagon in MRK. 
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A Dispersion relation, measure and pentagon transitions for gluons 


Here we review the building blocks of the Wilson loop OPE, coming from all gluonic exci¬ 
tations of the corresponding flux tube [29], in slightly adjusted notation. Although these 
formulas hold at finite coupling, they are particularly suited for the weak coupling expan¬ 
sion, as we will comment below, and display to leading order at the end of this section. 

The energy and the momentum of the a-th gluon bound state, a E Z* with a > 0 
corresponding to positive and a < 0 to negative helicity gluons, are given by 


E a {u) = |a| + Ag [Q • M • K(a,u)\ 1 , p a (u) = 2u - Ag [Q • M • n(a, u)] 1 , (A.l) 

where Q is a matrix with elements Qij = Sij(— l)* + 1 i, M is related to another matrix K , 


(1 + K)~ 1 = Y J (-K) n , 


Kij = 2.7 (— 


71=0 


oo 

2j(—i)A*+!) f dt Jj(2gt)Jj{‘2gt ) , ( A _ 2 ) 

J t e 1 


Ji is the i-th Bessel function of the first kind, and k, k are vectors with elements 

r dt Jj(2gt)(Jo(2gt) — cos(ut) [e*/ 2 ] ^ ^ ^ +1 ) 


e t — 1 


K 


;(a, u) = / -(-1) J+1 


dt . j + i Jj{2gt) sin(uf) [e t//2 ] ( 1J |a|+1 


(A.3) 


e t — 1 


The subscripts in (A.l) denote that we are only taking the first component of the vector 
inside the brackets. Moving on, the measure of the a-th gluon bound state, and its pentagon 
transition to a 6-th bound state are given by 

Ha(u) = F a (u)ef -3 

M a ’ b) (u,v)-ifi a ’ b) (u,v)+ft’ b) (n,v)-ft’ b) («,«) 


P a \b(u\v) = F atb (u,v)e l 


(A.4) 

(A.5) 
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where the f- a ' b ^ functions are expressed in terms of the previous quantities (A.2)-(A.3) 

f[ a,b ^ (u, v) = 2 k(a, u) ■ Q ■ M • n(b, v ), f^ b ) v) = 2 «(a, u) ■ Q ■ M • R(b, v) , 

ft’ b] (u, v) = 2 k(o, u) ■ Q • M • R(b, v ), /j a ’^ (it, v) = 2 /c(a, u) • Q • M • n(b, v), 

and in addition 

Jo" ( 2e ~ aV2 co S (ut)-M2gt)-l) 


(A.6) 


^a(lt) = 


< 7 2 (—1)“T(1 + | + m)r(l + | — iu)e Jo 


r(a)(x[+“]x[-“] - S 2 )V((® [+a] ) 2 - 5 2 )((x[-“]) 2 - 9 2 ) 


, (A.7) 


F a ,b(u, v ) = (.T[ +a ]y[ _b ] - 5r 2 )(x[ _a ]y[ +6 ] - g 2 )(a;[ +a V + ^ - g 2 )(®[ _a ]y[ _ ^ - g 2 ) x (A. 8) 

fOQ ^(J()(pg*)~l) / T n (o n f\ I 1 r — at/2 — iut r — bt/2+ivt\ 

(-l) fc T(+T +iu~iv)T{^ -iu + iv)e Jo ( o( 5 )+ ] 


F a _ 6 (n,i;) = 


S- 2 T(1 + | + m)T(l + | — iu)T(l + ^ — lit + in) 

rao dt(Jo(2gt)-l) ( Jg ( 2 g t ) +1 _ e - at / 2 -i u t_ e - bt /2+ivt) 

P, i(e 1) 


( 1 - 


■\+a]yl~b] 


^)(1- 


1 / [+b] 


)( 1 - 


xi+ a ]y[+ b ] 


)(! x [-a]yl-b ]) 


(A.9) 


r(l + + ra-w) 

T (1 + f + in)r(l + | - in) 


where for the last three formulas we have resticted a, b > 0, and under the same restriction 
the remaining cases may be obtained by F- a (u ) = F a (u), F_ a _b(u,v) = F a ^(u,v) and 
F_ atb (u,v) = F a _ b (u,v). Finally, 

x[ ±a ] = x(u ± ia/2) , = x(v ± ib/2) , x(u) = ^{u + ^/u 2 — (2g) 2 ). (A.10) 

Although the vectors (A.3) and matrices (A.3) are infinite-dimensional, and hence the 
matrix products (A.l) and (A.6) should really be thought of as infinite-dimensional sums 
at finite coupling, it is possible to show that at weak coupling Kjj starts at order 0 (g l+3 ). 
Thus if we wish to obtain the weak coupling expansion up to order 0(g 21 ), we may simply 
truncate all summations up to their first i. j = 1 , 2 ..., 21 — 1 terms. 

In order to see this, and also perform the expansion in practice, one starts with the 
Taylor series of the Bessel functions, 


J i( z ) = 


(-l) r 


—' k\T(i + 71 + 1) \2 

n=Q x 7 


z+2n 


and then computes integrals in t by means of 

t m e 


dt 


1 _- t ~ —Srnfi ) = (-ir + V m) (*) , m > 0, 


where is the polygamma function, 


fjm +1 

^ M W = ^+I logr( - ) 


(A.ll) 


(A.12) 


(A.13) 
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For example, the leading order expressions for all necessary ingredients of the Wilson loop 
OPE that we obtain in this manner, are for a, b > 0 


E a (u) — a + 2g 2 [ijj(l + ^ + iu) + '(/’(l + § — iu) — 2"0(1)] + 0(g 4 ) 
Pa(u) = 2u + 2ig 2 ('ip(% + iu) - ?/>(§ - iu)) + 0(g 4 ), 


A»a(«) = (-l)V 


P(| + w)P(§ ~ w) 

r(a)(«+f)(«-f) 


+ 0(<7 4 




(-1) 6 (| - iu){\ + iv)T{^ + iu- iv)T( 2±6 - m + iu) 


+ 0(2° 


(A.14) 

(A.15) 

(A.16) 

(A.17) 


g 2 r(| + m)r(| — iu)r(l + ^2 - iu + i y ) 

More details about the weak coupling expansion may be found in appendix E of [26], as 
well as appendix A.2 of [29]). 


B Remarks on the two-particle gluon bound state contributions 

Although we leave the resummation of all two-particle gluon bound states as an exciting 
open question for future work, in this appendix we shall employ the technology we have 
developed so far in order to compute the contribution of several such states in the collinear 
limit (part B.l). The main motivation here is to examine whether these states survive 
when we transition from collinear to multi-Regge kinematics, and complement the analysis 
of section 3 (part B.2) . We gather evidence up to 5 loops that this is not the case, thus 
allowing us to reconstruct the hexagon remainder function in the latter kinematics from 
the single-particle gluon bound states, resummed in subsection 2.3. 


B.l Evaluation of the integrals 

The two-particle gluon bound state OPE contribution corresponds to keeping only the 
N = 2 term in eq. (2.13), which we may rewrite as 

WWm ... f + °° [ + °° dudv p, a (u)p,b(v) 




V 2 [aM 

1+5, 


ab 


^ 2 la,b] = 


' — oo «/ —oo 


(2vr) 2 P a \ b {u\v)P b \ a (v\u) ’ 


(B.l) 


a= 1 6=1 

with a, b the number of gluons contained in each of the two bound states (also equal to their 
twist and helicity), 5 a b a Kronecker delta function, and all ingredients of the OPE integral 
on the right-hand side contained in appendix A. With the help of (A.14)-(A.17), at 4 loops 
the integral becomes 

■ dudv r (| - iu) 2 r (| + iu) 2 r (| - iv) 2 r (§ + iv) 2 


(B.2) 


yy( 4 ) = 

2 Kb] j ( 2vr )2 r( a )r(6)r (§ +§ + iu - iv) r (§ + \ - iu + iv) " 

(-l) a ~ b [+ - i{u - tp] [gr + i( U - + (a+6)(r+^)+2i«r(«+ v) 

(u-f)2(u+f)2(u-f)2(u+f)2 

where we have replaced four of the gamma functions coming from (A.17) with the help of 
the identity 

r (l + ^ - iu + iv) T (l + ^ - iv + iu) ,^ a -b[a-b -t m ra-6 ■/ m 

T(*=± + iu-iv)r(^+iv-iu) K 12 v n L 2 v yj ’ 

(B.3) 


21 











which may be derived from (2.17). For specific values of a, 6, and depending on whether 
a/2,6/2 and (a + 6)/2 are integer or half-integer, we may similarly replace the remaining 
gamma functions by virtue of 


T(n + x)T(n — x) = TT (m + x)(m — x) — 7 —r , 

- LJ - sm(7rx) 

m= 1 

71—1 ^ 

r(i + n + x)r(i +n — x)= TT (i + m + x)(b + m — x ) - 7 — T , 

- LJ - z z cos ( ttx ) 

m =0 v ’ 


(B.4) 


for n integer. Finally, we may factorise all u— and v —dependence coming from the trigono¬ 
metric/hyperbolic functions, e.g. 


sinh( 7 ru — irv) = sinh(7rit) cosh(7rn) — cosh(7ru) sinh(7rn). 


(B.5) 


From these considerations, it is clear that the only inseparable dependence on u — v 
will come from the products in (B.4) appearing in the denominator of (B.l). Given that 
at higher loops the integrand (B.l) is dressed by sums of products of polygamma functions 
with arguments 1 + |±«, |±it, l + |=tu and |± v, this statement holds true independently 
of the loop order. 

As was first observed in [31], for the case a = 6 = 1 these (it — independent terms in 
the denominator cancel out, and therefore VV^jj ^ always reduces to a sum of two factorised 
1-fold integrals, which may be evaluated along the lines of [30]. Here, we find that exactly 
the same phenomenon occurs also for a = 2 and 6=1. 

Although for a + 6 > 4 we start dealing with genuine 2-fold integrals, these can still be 
done by summing over residues, at least at 4 loops. The main subtlety is that we need to 
treat separately not only the poles at it = ia /2 and v = ib/2 , but also the ones where 

u - v = *(-2±$ + l),i(-“±$ + 2),...,t(^-l) , u- v ^±i^. (B.6) 

These special poles will lead to simple sums, whereas the remaining ones give rise to double 
sums, which may be brought to a form similar to the one we encountered in 2.3. Sparing 
the reader the rest of the details, all in all we arrive at expressions for ^ and VV,^ ^ 

up to l = 5 loops, as well W,^ ^ and 2 ] • which are included in the accompanying hie 

boundstatesN2.m. 

Before moving on to examine the contribution of these states in MRK, let us end this 
section by noting a peculiar property of the two twist-4 contributions: Each one separately 
gives rise to polylogarithms with symbol letters 5,1 + S 2 but also 1 — S' 2 , namely HPLs with 
both positive and negative weights. At first this may seem to contradict the expectation 
discussed in section 2.4, that perturbatively Rq, and hence also W, are described by a 
particular 9-letter alphabet, since the latter only reduces to the letters S and 1 + S 2 in the 
collinear limit. Quite remarkably however, in their sum ^ + W.^ 9 j/ 2, also taking into 
account the symmetry factor in B.l, all terms containing the letter 1 — S 2 cancel out, and 
consistency is restored. This perhaps suggests that it is only the combination of states with 
the same charges and particle number N that has a physical significance in the OPE. 
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B.2 Contribution to MRK 

In section 3.1, we reviewed the analytic continuation which is necessary for obtaining a 
nontrivial result in MRK, and specialised its form in a regime where it overlaps with the 
double scaling limit, where our resummed single-particle gluon contribution lives. For the 
individual two-particle gluon bound states that we computed in the previous part of this 
appendix, we will need a similar specialisation to a region overlapping with the collinear 
limit. This was first considered in [62], based on the initial Wilson loop OPE approach 
[23]. Following the refinement of the latter approach, this procedure has been revisited and 
extended in [36], which we now briefly review. 

In the latter reference, the usual analytic continuation in the cross-ratios was translated 
to T,a,<p variables (2.3) parameterising the near-collinear limit expansion. In more detail, 
if we denote S = e a and T = e~ T , a particular path 

C : (e iX S, T, e lx cos 4> + z<S ,_1 (T + T _1 ) sin^), x G [0,7r], (B.7) 

connecting the initial point ( S , T, cos 4>) in the Euclidean sheet for x = 0, to the same point 
(-S, T, — cos 4>) in the 2 —> 4 Mandelstam sheet for % = tv. 

In [36] it was also observed that at two loops, the analytic continuation for S > 1 and 
the collinear limit expansion commute. In other words, first analytically continuing 
and then expanding the near collinear limit yields the same result as first expanding, and 
then analytically continuing term by term. This behaviour was then conjectured to hold 
to all loops, and in fact it can be proven for all polylogarithms on Ato,6- In other words if 
one assumes that Rq is described by hexagon functions, then the commuting of the two 
procedures follows. 

Taking the aforementioned property as granted, it is then easy to obtain the collinear 
limit expansion of the analytically continued R$. For all (sums of) positive helicity gluon 
states with given particle number and twist a that we have considered (as discussed in the 
last paragraph of section B.l), we find that their kinematical dependence is always a linear 
combination of terms of the form 

(logT) m f — J S 2k (logSyH mu ..., mr (-l/S 2 ) , k = 0,1... a , m > 0 . (B.8) 

For the rational part in and S, clearly the value at the endpoint of the analytic continu¬ 
ation does not depend on the path, so that we can simply set e 1 ^ —> —e* 1 ^, S —> —S, and in 
fact all extra minus signs always cancel out. The HPLs with positive weights also remain 
unchanged, since they only have branch cuts for their argument between 1 and infinity, 
and the path (B.7) never crosses them for S > 0. Therefore analytically continuing our 
expressions amounts to the almost trivial replacement 

log S —> log S + in , or a —> a + in . (B.9) 

The final step is to pass from the analytically continued collinear to multi-Regge kinematics. 
As was first noticed in [16], and can be readily verified from (2.3) and (3.1), in the variables 
we’re currently using MRK corresponds to T —> 0, S —>• 0 with r = T/S fixed. The ratio 
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r has already appeared in (3.3), and it is in fact via this procedure that the relation with 
the w, w* variables of the MRK appearing on the left-hand side of the latter formula can 
be established. 

We thus arrive at the following connection between the two kinematics regions: Start¬ 
ing with the analytically continued collinear limit expansion (around T = 0), expanding 
additionally around S = 0 and replacing T = rS lands us on the multi-Regge limit expan¬ 
sion (around 5 = 0), where we are additionally expanding around r = 0. It is specifically 
the non-power suppressed 0(5°) term that has principally been the focus of BFKL analysis, 
and also of our paper here. 

For that reason, we will take the strict 5 —> 0 limit of our expressions (up to large 
logarithms), under which we see that, very interestingly, all but the k = 0 part of all 
terms (B.8) already drops out. Further examining the k = 0 term for W^j 1 j, ^ up to 

l = 5, and also ji + W,^ 2]/2, we see that it vanishes in all cases. We deem this as 
very compelling evidence that 2-particle gluon bound states do not contribute in MRK in 
general, since this behaviour remain unchanged across loop orders, and also for qualitatively 
quite different types of integrals. 

C Multiple Polylogarithms 

In this appendix we review the definitions and several properties of multiple polylogarithms 
(MPLs), which will be useful throughout the main text. More detailed expositions may be 
found for example in [63-65]. 

Multiple polylogarithms may be defined as nested sums 11 

“ x * 1 x ik 

Lim ll ... 1 m fc (®l) ■ ■ • , Xk) = / ] • • • vjjf J (C-l) 

*l>*2>...>*fe>0 k 

in the region 


\x\X 2 ■ ■ ■ Xj\ < 1 Vj = l,...k, with (mi, x\) ^ (1,1), (C.2) 

so that the above series converges. Comparing (C.l) with the definition of Z-sums (2.21), 
it is evident that MPLs are special cases of the latter with the outer summation index 
set to infinity, Li {xi, • • • , Xj ) = Z(oo; m \,..., my,x i,..., Xj). Thus similarly to the 
discussion of Z-surns around eq. (2.22), rearranging the summation ranges of a product of 
MPLs gives rise to quasi-shuffle relations equating the product to a linear combination of 
MPLs. For example, 

Limi (xi )Li m2 ( 3 ^ 2 ) — Li mi)m2 {x \, X 2 ) T Li m2 mi (x 2 , x\) -f- ( X 1 X 2 ) • (C.3) 

The nested sum definition (C.l) is more suited for the evaluation of Mellin-Bernes type 
integrals by residues, as we did in section 2. However there also exists an alternative, global 

11 Note that a different convention with the order of the summation indices reversed, i.e. ik > ■ ■ ■ > ii, 
also exists in the literature. 
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definition of MPLs in terms of (regularised) iterated integrals, as 


G(oi, ...,a n ;z)= < 




if a\ = ... = a n = 0 


,io tp^T G '( a 2 ) • ■ ■ ,a n ;ti ), G(;z) = 1, otherwise. 

Unfolding the generic case of the second line for 7 ^ 0, we may equivalently write it as 

dt\ [ ll dt2 f tn ~ i dt n 


G{ai, ...,a n ;z) = 


tl ~ Or J o f2 — 


t n CL 11 


(C.4) 

as 

(C.5) 


a = (ai,..., a n ) is usually denoted as the singularity vector, and its length n as the weight 
or transcendentality of the MPL. 

Although seemingly different, the sum and integral definitions are equivalent in the 
region (C.2), and in particular 


Limi,...,m fe (®ij — ( 1) ' G ( 0, 0 , ^ , 0, , 0 , x x ' 0-"; 0 j ^ ^ ■ 

mi —1 m 2 —1 m/g —1 


1). (C. 6 ) 


A useful identity which follows immediately from (C.5), is that for 7 ^ 0 G-functions are 
invariant under an arbitrary rescaling x G C* of all the arguments, namely 


G(ai,..., a^; z) = G(xai, ...,xa k ;xz). 


(C.7) 


Finally, the integral dehnition of MPLs reveals the existence of further shiijfle algebra 
relations between products of G-functions with the same rightmost argument, such as 


G{a-z)G(b;z) = 


dt\ 

h- a 


dt2 
t'2 - b 
dt\ dt2 

lo<t 2 <ti<z ti — a t 2 — b 
= G(a, 6 ; z) + G(5, a; z). 


+ 


dt 1 dt2 

0 <t!<t 2 <z ti — a t 2 — b 


(C. 8 ) 


More generally, 

G{a;z)G(b;z) = ^ G(c;z) (C.9) 

cEa III b 

where the shuffle product a HI b is defined as the set of all permutations of the elements of 
aU b, that preserve the ordering among the elements a, £ a, and among the elements £ b. 

Interesting special cases of MPLs include the harmonic polylogarithms [46], for singu¬ 
larity vector entries dj £ {0,±1}, and the two-dimensional harmonic polylogarithms [35], 
for ai £ { 0 , 1 , — w, 1 — w}, a subset of which we have encountered in this paper 12 . 

12 In the notation of [35], 2dHPLs are also defined with an extra minus sign when the sum of all the ai 
equal to 1 or 1 — 2 is odd. 


25 














References 


[1] L. F. Alday and J. M. Maldacena, Gluon scattering amplitudes at strong coupling, JHEP 
0706 (2007) 064, [arXiv: 0705. 0303], 

[2] J. Drummond, G. Korchemsky, and E. Sokatchev, Conformal properties of four-gluon planar 
amplitudes and Wilson loops, Nucl.Phys. B795 (2008) 385-408, [arXiv: 0707. 0243]. 

[3] A. Brandhuber, P. Heslop, and G. Travaglini, MHV amplitudes in N = 4 super Yang-Mills 
and Wilson loops, Nucl.Phys. B794 (2008) 231-243, [arXiv: 0707. 1153]. 

[4] J. Drummond, J. Henn, G. Korchemsky, and E. Sokatchev, On planar gluon 
amplitudes/Wilson loops duality, Nucl.Phys. B795 (2008) 52-68, [arXiv: 0709. 2368]. 

[5] Z. Bern, L. Dixon, D. Kosower, R. Roiban, M. Spradlin, et al., The Two-Loop Six-Gluon 
MHV Amplitude in Maximally Supersymmetric Yang-Mills Theory, Phys.Rev. D78 (2008) 
045007, [arXiv: 0803. 1465], 

[6] J. Drummond, J. Henn, G. Korchemsky, and E. Sokatchev, Hexagon Wilson loop = six-gluon 
MHV amplitude, Nucl.Phys. B815 (2009) 142-173, [arXiv: 0803. 1466]. 

[7] N. Berkovits and J. Maldacena, Fermionic T-Duality, Dual Superconformal Symmetry, and 
the Amplitude/Wilson Loop Connection, JHEP 0809 (2008) 062, [arXiv:0807.3196[. 

[8] L. Mason and D. Skinner, The Complete Planar S-matrix of A/” = 4 SYM as a Wilson Loop 
in Twistor Space, JHEP 1012 (2010) 018, [arXiv: 1009. 2225]. 

[9] S. Caron-Huot, Notes on the scattering amplitude / Wilson loop duality, JHEP 1107 (2011) 
058, [arXiv: 1010.1167], 

[10] J. Drummond, J. Henn, G. Korchemsky, and E. Sokatchev, Conformal Ward identities for 
Wilson loops and a test of the duality with gluon amplitudes, Nucl.Phys. B826 (2010) 
337-364, [arXiv: 0712. 1223], 

[11] V. Del Duca, C. Duhr, and V. A. Smirnov, The Two-Loop Hexagon Wilson Loop mAl = 4 
SYM, JHEP 1005 (2010) 084, [arXiv: 1003. 1702], 

[12] A. B. Goncharov, M. Spradlin, C. Vergu, and A. Volovich, Classical Polylogarithms for 
Amplitudes and Wilson Loops, Phys.Rev.Lett. 105 (2010) 151605, [arXiv: 1006. 5703]. 

[13] S. Caron-Huot, Superconformal symmetry and two-loop amplitudes in planar N = 4 super 
Yang-Mills, JHEP 1112 (2011) 066, [arXiv: 1105.5606]. 

[14] L. J. Dixon, J. M. Drummond, and J. M. Henn, Bootstrapping the three-loop hexagon, JHEP 
1111 (2011) 023, [arXiv: 1108.4461], 

[15] L. J. Dixon, J. M. Drummond, and J. M. Henn, Analytic result for the two-loop six-point 
NMHV amplitude in AT = 4 super Yang-Mills theory, JHEP 1201 (2012) 024, 

[arXiv: 1111. 1704], 

[16] L. J. Dixon, J. M. Drummond, M. von Hippel, and J. Pennington, Hexagon functions and 
the three-loop remainder function, JHEP 1312 (2013) 049, [arXiv: 1308.2276]. 

[17] L. J. Dixon, J. M. Drummond, C. Duhr, and J. Pennington, The four-loop remainder 
function and multi-Regge behavior at NNLLA in planar JN = 4 super-Yang-Mills theory, 
JHEP 1406 (2014) 116, [arXiv: 1402. 3300], 

[18] L. J. Dixon and M. von Hippel, Bootstrapping an NMHV amplitude through three loops, 
JHEP 1410 (2014) 65, [arXiv: 1408.1505], 


26 


[19] J. M. Drummond, G. Papathanasiou, and M. Spradlin, A Symbol of Uniqueness: The Cluster 
Bootstrap for the 3-Loop MHV Heptagon, JHEP 1503 (2015) 072, [arXiv: 1412. 3763]. 

[20] J. Golden, A. B. Goncharov, M. Spradlin, C. Vergu, and A. Volovich, Motivic Amplitudes 
and Cluster Coordinates, JHEP 1401 (2014) 091, [arXiv: 1305.1617]. 

[21] S. Fomin and A. Zelevinsky, Cluster algebras I: Foundations, J. Am. Math. Soc. 15 (2002), 
no. 2 497-529, [math/0104151]. 

[22] S. Fomin and A. Zelevinsky, Cluster algebras II: Finite type classification, Invent. Math. 154 
(2003), no. 1 63-121, [math/0208229]. 

[23] L. F. Alday, D. Gaiotto, J. Maldacena, A. Sever, and P. Vieira, An Operator Product 
Expansion for Polygonal null Wilson Loops, JHEP 1104 (2011) 088, [arXiv: 1006.2788]. 

[24] D. Gaiotto, J. Maldacena, A. Sever, and P. Vieira, Pulling the straps of polygons, JHEP 
1112 (2011) 011, [arXiv: 1102.0062], 

[25] B. Basso, A. Sever, and P. Vieira, Spacetime and Flux Tube S-Matrices at Finite Coupling 
for J\T = 4 Supersymmetric Yang-Mills Theory, Phys.Rev.Lett. Ill (2013), no. 9 091602, 
[arXiv: 1303.1396], 

[26] B. Basso, A. Sever, and P. Vieira, Space-time S-matrix and Flux tube S-matrix II. Extracting 
and Matching Data, JHEP 1401 (2014) 008, [arXiv: 1306.2058]. 

[27] B. Basso, A. Sever, and P. Vieira, Space-time S-matrix and Flux-tube S-matrix III. The 
two-particle contributions, JHEP 1408 (2014) 085, [arXiv : 1402. 3307]. 

[28] B. Basso, Exciting the GKP string at any coupling, Nucl.Phys. B857 (2012) 254-334, 

[arXiv: 1010.5237], 

[29] B. Basso, A. Sever, and P. Vieira, Space-time S-matrix and Flux-tube S-matrix IV. Gluons 
and Fusion, JHEP 1409 (2014) 149, [arXiv: 1407. 1736]. 

[30] G. Papathanasiou, Hexagon Wilson Loop OPE and Harmonic Polylogarithms, JHEP 1311 
(2013) 150, [arXiv: 1310.5735], 

[31] G. Papathanasiou, Evaluating the six-point remainder function near the collinear limit, 

Int.J.Mod.Phys. A29 (2014), no. 27 1450154, [arXiv: 1406. 1123]. 

[32] S. Moch, P. Uwer, and S. Weinzierl, Nested sums, expansion of transcendental functions and 
multiscale multiloop integrals, J.Math.Phys. 43 (2002) 3363-3386, [hep-ph/0110083]. 

[33] S. Weinzierl, Symbolic expansion of transcendental functions, Comput.Phys.Commun. 145 
(2002) 357-370, [math-ph/0201011], 

[34] F. C. Brown, Multiple zeta values and periods of moduli spaces 91to. n (R), Annales Sci.Ecole 
Norm.Sup. 42 (2009) 371, [math/0606419]. 

[35] T. Gehrmann and E. Remiddi, Two loop master integrals for gamma* —> 3 jets: The 
Planar topologies, Nucl.Phys. B601 (2001) 248-286, [hep-ph/0008287]. 

[36] Y. Hatsuda, Wilson loop OPE, analytic continuation and multi-Regge limit, JHEP 1410 
(2014) 38, [arXiv: 1404.6506], 

[37] L. J. Dixon, C. Duhr, and J. Pennington, Single-valued harmonic polylogarithms and the 
multi-Regge limit, JHEP 1210 (2012) 074, [arXiv: 1207.0186]. 

[38] B. Basso, S. Caron-Huot, and A. Sever, Adjoint BFKL at finite coupling: a short-cut from 
the collinear limit, arXiv: 1407.3766. 


27 


[39] A. Belitsky, Nonsinglet pentagons and NHMV amplitudes, arXiv: 1407.2853. 

[40] A. Belitsky, Fermionic pentagons and NMHV hexagon, arXiv: 1410.2534. 

[41] B. Basso, J. Caetano, L. Cordova, A. Sever, and P. Vieira, OPE for all Helicity Amplitudes, 
arXiv:1412.1132. 

[42] A. V. Belitsky, On factorization of multiparticle pentagons, Nucl. Phys. B897 (2015) 
346-373, [arXiv: 1501.0686], 

[43] C. Bauer, A. Frink, and R. Kreckel, Introduction to the ginac framework for symbolic 
computation within the C++ programming language, J.Symb.Comput. 33 (2002) 1-12, 
[cs/0004015], 

[44] S. Moch and P. Uwer, -xsummer- transcendental functions and symbolic summation in form, 
Computer Physics Communications 174 (2006), no. 9 759 - 770. 

[45] J. A. M. Vermaseren, New features of FORM, math-ph/0010025. 

[46] E. Remiddi and J. Vermaseren, Harmonic polylogarithms, Int.J.Mod.Phys. A15 (2000) 
725-754, [hep-ph/9905237], 

[47] L. Dixon. Private communication, August 2014. 

[48] J. Bartels, L. Lipatov, and A. Sabio Vera, BFKL Pomeron, Reggeized gluons and 
Bern-Dixon-Smirnov amplitudes, Phys.Rev. D80 (2009) 045002, [arXiv: 0802. 2065]. 

[49] J. Bartels, L. Lipatov, and A. Sabio Vera, AT = 4 supersymmetric Yang Mills scattering 
amplitudes at high energies: The Regge cut contribution, Eur.Phys.J. C65 (2010) 587-605, 
[arXiv: 0807. 0894], 

[50] L. Lipatov and A. Prygarin, Mandelstam cuts and light-like Wilson loops in AT = 4 SUSY, 
Phys.Rev. D83 (2011) 045020, [arXiv: 1008. 1016], 

[51] L. Lipatov and A. Prygarin, BFKL approach and six-particle MHV amplitude in AT = 4 
super Yang-Mills, Phys.Rev. D83 (2011) 125001, [arXiv: 1011. 2673]. 

[52] J. Bartels, L. Lipatov, and A. Prygarin, MHV Amplitude for 3—>3 Gluon Scattering in Regge 
Limit, Phys.Lett. B705 (2011) 507-512, [arXiv: 1012. 3178]. 

[53] V. Fadin and L. Lipatov, BFKL equation for the adjoint representation of the gauge group in 
the next-to-leading approximation at AT = 4 SUSY, Phys.Lett. B706 (2012) 470-476, 

[arXiv: 1111.0782], 

[54] A. Prygarin, M. Spradlin, C. Vergu, and A. Volovich, All Two-Loop MHV Amplitudes in 
Multi-Regge Kinematics From Applied Symbology, Phys.Rev. D85 (2012) 085019, 

[arXiv: 1112.6365], 

[55] L. Lipatov, A. Prygarin, and H. J. Schnitzer, The Multi-Regge limit of NMHV Amplitudes in 
AT = 4 SYM Theory, JHEP 1301 (2013) 068, [arXiv: 1205.0186], 

[56] J. Bartels, A. Kormilitzin, L. Lipatov, and A. Prygarin, BFKL approach and 2 —» 5 
maximally helicity violating amplitude in AT = 4 super-Yang-Mills theory, Phys.Rev. D86 
(2012) 065026, [arXiv: 1112. 6366], 

[57] J. Bartels, A. Kormilitzin, and L. Lipatov, Analytic structure of the n = 7 scattering 
amplitude in Af = 4 SYM theory at multi-Regge kinematics: Conformal Regge pole 
contribution, Phys.Rev. D89 (2014) 065002, [arXiv: 1311 .2061]. 


28 


[58] J. Bartels, A. Kormilitzin, and L. N. Lipatov, Analytic structure of the n = 7 scattering 
amplitude in Af = 4 SYM theory in multi-Regge kinematics: Conformal Regge cut 
contribution , arXiv: 1411.2294. 

[59] D. Maitre, HPL, a mathematica implementation of the harmonic polylogarithms, 
Comput.Phys.Commun. 174 (2006) 222-240, [hep-ph/0507152]. 

[60] F. Brown, Poly log arithmes multiples uniformes en une variable, C.R. Acad. Sci. Paris Ser. I 
338 (2004) 527-532. 

[61] J. Drummond, C. Duhr, B. Eden, P. Heslop, J. Pennington, et al., Leading singularities and 
off-shell conformal integrals, JHEP 1308 (2013) 133, [arXiv : 1303. 6909]. 

[62] J. Bartels, L. Lipatov, and A. Prygarin, Collinear and Regge behavior of 2 -> f MHV 
amplitude in N = f super Yang-Mills theory , arXiv: 1104.4709. 

[63] S. Weinzierl, Introduction to Feynman Integrals, in 6th Summer School on Geometric and 
Topological Methods for Quantum Field Theory Villa de Leyva, Colombia, July 6-23, 2009, 
2010. arXiv:1005.1855. 

[64] C. Duhr, Mathematical aspects of scattering amplitudes, arXiv: 1411.7538. 

[65] E. Panzer, Feynman integrals and hyperlogarithms, arXiv: 1506.0724. 


29 


